Variant-to-function analysis of the childhood obesity chr12q13 locus implicates rs7132908 as a causal variant within the 3' UTR of FAIM2.

The ch12q13 locus is among the most significant childhood obesity loci identified in genome-wide association studies. This locus resides in a non-coding region within FAIM2; thus, the underlying causal variant(s) presumably influence disease susceptibility via cis-regulation. We implicated rs7132908 as a putative causal variant by leveraging our in-house 3D genomic data and public domain datasets. Using a luciferase reporter assay, we observed allele-specific cis-regulatory activity of the immediate region harboring rs7132908. We generated isogenic human embryonic stem cell lines homozygous for either rs7132908 allele to assess changes in gene expression and chromatin accessibility throughout a differentiation to hypothalamic neurons, a key cell type known to regulate feeding behavior. The rs7132908 obesity risk allele influenced expression of FAIM2 and other genes and decreased the proportion of neurons produced by differentiation. We have functionally validated rs7132908 as a causal obesity variant that temporally regulates nearby effector genes and influences neurodevelopment and survival.

The association between neighborhood deprivation and DNA methylation in an autopsy cohort.

Previous research has found that living in a disadvantaged neighborhood is associated with poor health outcomes. Living in disadvantaged neighborhoods may alter inflammation and immune response in the body, which could be reflected in epigenetic mechanisms such as DNA methylation (DNAm). We used robust linear regression models to conduct an epigenome-wide association study examining the association between neighborhood deprivation (Area Deprivation Index; ADI), and DNAm in brain tissue from 159 donors enrolled in the Emory Goizueta Alzheimer's Disease Research Center (Georgia, USA). We found one CpG site (cg26514961, gene PLXNC1) significantly associated with ADI after controlling for covariates and multiple testing (p-value=5.0e-8). Effect modification by APOE ε4 was statistically significant for the top ten CpG sites from the EWAS of ADI, indicating that the observed associations between ADI and DNAm were mainly driven by donors who carried at least one APOE ε4 allele. Four of the top ten CpG sites showed a significant concordance between brain tissue and tissues that are easily accessible in living individuals (blood, buccal cells, saliva), including DNAm in cg26514961 (PLXNC1). Our study identified one CpG site (cg26514961, PLXNC1 gene) that was significantly associated with neighborhood deprivation in brain tissue. PLXNC1 is related to immune response, which may be one biological pathway how neighborhood conditions affect health. The concordance between brain and other tissues for our top CpG sites could make them potential candidates for biomarkers in living individuals.

Epigenetic associations in HPA axis genes related to bronchopulmonary dysplasia and antenatal steroids.

BACKGROUND: Bronchopulmonary dysplasia (BPD), a common morbidity among very preterm infants, is associated with chronic disease and neurodevelopmental impairments. A hypothesized mechanism for these outcomes lies in altered glucocorticoid (GC) activity. We hypothesized that BPD and its treatments may result in epigenetic differences in the hypothalamic-pituitary-adrenal (HPA) axis, which is modulated by GC, and could be ascertained using an established GC risk score and DNA methylation (DNAm) of HPA axis genes. METHODS: DNAm was quantified from buccal tissue (ECHO-NOVI) and from neonatal blood spots (ELGAN ECHO) via the EPIC microarray. Prenatal maternal characteristics, pregnancy complication, and neonatal medical complication data were collected from medical record review and maternal interviews. RESULTS: The GC score was not associated with steroid exposure or BPD. However, six HPA genes involved in stress response regulation demonstrated differential methylation with antenatal steroid exposure; two CpGs within FKBP5 and POMC were differentially methylated with BPD severity. These findings were sex-specific in both cohorts; males had greater magnitude of differential methylation within these genes. CONCLUSIONS: These findings suggest that BPD severity and antenatal steroids are associated with DNAm at some HPA genes in very preterm infants and the effects appear to be sex-, tissue-, and age-specific. IMPACT: This study addresses bronchopulmonary dysplasia (BPD), an important health outcome among preterm neonates, and interrogates a commonly studied pathway, the hypothalamic-pituitary-adrenal (HPA) axis. The combination of BPD, the HPA axis, and epigenetic markers has not been previously reported. In this study, we found that BPD itself was not associated with epigenetic responses in the HPA axis in infants born very preterm; however, antenatal treatment with steroids was associated with epigenetic responses.

Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits.

The hypothalamus regulates metabolic homeostasis by influencing behavior and endocrine systems. Given its role governing key traits, such as body weight and reproductive timing, understanding the genetic regulation of hypothalamic development and function could yield insights into disease pathogenesis. However, given its inaccessibility, studying human hypothalamic gene regulation has proven challenging. To address this gap, we generate a high-resolution chromatin architecture atlas of an established embryonic stem cell derived hypothalamic-like neuron model across three stages of in vitro differentiation. We profile accessible chromatin and identify physical contacts between gene promoters and putative cis-regulatory elements to characterize global regulatory landscape changes during hypothalamic differentiation. Next, we integrate these data with GWAS loci for various complex traits, identifying multiple candidate effector genes. Our results reveal common target genes for these traits, potentially affecting core developmental pathways. Our atlas will enable future efforts to determine hypothalamic mechanisms influencing disease susceptibility.

CRISPR-Cas9-Mediated Genome Editing Confirms EPDR1 as an Effector Gene at the BMD GWAS-Implicated 'STARD3NL' Locus.

Genome-wide-association studies (GWASs) have discovered genetic signals robustly associated with BMD, but typically not the precise localization of effector genes. By intersecting genome-wide promoter-focused Capture C and assay for transposase-accessible chromatin using sequencing (ATAC-seq) data generated in human mesenchymal progenitor cell (hMSC)-derived osteoblasts, consistent contacts were previously predicted between the EPDR1 promoter and multiple BMD-associated candidate causal variants at the 'STARD3NL' locus. RNAi knockdown of EPDR1 expression in hMSC-derived osteoblasts was shown to lead to inhibition of osteoblastogenesis. To fully characterize the physical connection between these putative noncoding causal variants at this locus and the EPDR1 gene, clustered regularly interspaced short-palindromic repeat Cas9 endonuclease (CRISPR-Cas9) genome editing was conducted in hFOB1.19 cells across the single open-chromatin region harboring candidates for the underlying causal variant, rs1524068, rs6975644, and rs940347, all in close proximity to each other. RT-qPCR and immunoblotting revealed dramatic and consistent downregulation of EPDR1 specifically in the edited differentiated osteoblast cells. Consistent with EPDR1 expression changes, alkaline phosphatase staining was also markedly reduced in the edited differentiated cells. Collectively, CRISPR-Cas9 genome editing in the hFOB1.19 cell model supports previous observations, where this regulatory region harboring GWAS-implicated variation operates through direct long-distance physical contact, further implicating a key role for EPDR1 in osteoblastogenesis and BMD determination. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Biological constraints on GWAS SNPs at suggestive significance thresholds reveal additional BMI loci.

To uncover novel significant association signals (p<5×10-8), genome-wide association studies (GWAS) requires increasingly larger sample sizes to overcome statistical correction for multiple testing. As an alternative, we aimed to identify associations among suggestive signals (5 × 10-8≤p<5×10-4) in increasingly powered GWAS efforts using chromatin accessibility and direct contact with gene promoters as biological constraints. We conducted retrospective analyses of three GIANT BMI GWAS efforts using ATAC-seq and promoter-focused Capture C data from human adipocytes and embryonic stem cell (ESC)-derived hypothalamic-like neurons. This approach, with its extremely low false-positive rate, identified 15 loci at p<5×10-5 in the 2010 GWAS, of which 13 achieved genome-wide significance by 2018, including at NAV1, MTIF3, and ADCY3. Eighty percent of constrained 2015 loci achieved genome-wide significance in 2018. We observed similar results in waist-to-hip ratio analyses. In conclusion, biological constraints on sub-significant GWAS signals can reveal potentially true-positive loci for further investigation in existing data sets without increasing sample size.

Variant-to-Gene-Mapping Analyses Reveal a Role for the Hypothalamus in Genetic Susceptibility to Inflammatory Bowel Disease.

BACKGROUND & AIMS: Inflammatory bowel disease (IBD) is a polygenic disorder characterized principally by dysregulated inflammation impacting the gastrointestinal tract. However, there also is increasing evidence for a clinical association with stress and depression. Given the role of the hypothalamus in stress responses and in the pathogenesis of depression, useful insights could be gleaned from understanding its genetic role in IBD. METHODS: We conducted genetic correlation analyses on publicly available genome-wide association study summary statistics for depression and IBD traits to identify genetic commonalities. We used partitioned linkage disequilibrium score regression, leveraging our ATAC sequencing and promoter-focused Capture C data, to measure enrichment of IBD single-nucleotide polymorphisms within promoter-interacting open chromatin regions of human embryonic stem cell-derived hypothalamic-like neurons (HNs). Using the same data sets, we performed variant-to-gene mapping to implicate putative IBD effector genes in HNs. To contrast these results, we similarly analyzed 3-dimensional genomic data generated in epithelium-derived colonoids from rectal biopsy specimens from donors without pathologic disease noted at the time of colonoscopy. Finally, we conducted enrichment pathway analyses on the implicated genes to identify putative IBD dysfunctional pathways. RESULTS: We found significant genetic correlations (rg) of 0.122 with an adjusted P (Padj) = 1.4 × 10-4 for IBD: rg = 0.122; Padj = 2.5 × 10-3 for ulcerative colitis and genetic correlation (rg) = 0.094; Padj = 2.5 × 10-3 for Crohn's disease, and significant approximately 4-fold (P = .005) and approximately 7-fold (P = .03) enrichment of IBD single-nucleotide polymorphisms in HNs and colonoids, respectively. We implicated 25 associated genes in HNs, among which CREM, CNTF, and RHOA encode key regulators of stress. Seven genes also additionally were implicated in the colonoids. We observed an overall enrichment for immune and hormonal signaling pathways, and a colonoid-specific enrichment for microbiota-relevant terms. CONCLUSIONS: Our results suggest that the hypothalamus warrants further study in the context of IBD pathogenesis.

Mapping effector genes at lupus GWAS loci using promoter Capture-C in follicular helper T cells.

Systemic lupus erythematosus (SLE) is mediated by autoreactive antibodies that damage multiple tissues. Genome-wide association studies (GWAS) link >60 loci with SLE risk, but the causal variants and effector genes are largely unknown. We generated high-resolution spatial maps of SLE variant accessibility and gene connectivity in human follicular helper T cells (TFH), a cell type required for anti-nuclear antibodies characteristic of SLE. Of the ~400 potential regulatory variants identified, 90% exhibit spatial proximity to genes distant in the 1D genome sequence, including variants that loop to regulate the canonical TFH genes BCL6 and CXCR5 as confirmed by genome editing. SLE 'variant-to-gene' maps also implicate genes with no known role in TFH/SLE disease biology, including the kinases HIPK1 and MINK1. Targeting these kinases in TFH inhibits production of IL-21, a cytokine crucial for class-switched B cell antibodies. These studies offer mechanistic insight into the SLE-associated regulatory architecture of the human genome.

Type 1 diabetes in Africa: an immunogenetic study in the Amhara of North-West Ethiopia.

AIMS/HYPOTHESIS: We aimed to characterise the immunogenic background of insulin-dependent diabetes in a resource-poor rural African community. The study was initiated because reports of low autoantibody prevalence and phenotypic differences from European-origin cases with type 1 diabetes have raised doubts as to the role of autoimmunity in this and similar populations. METHODS: A study of consecutive, unselected cases of recently diagnosed, insulin-dependent diabetes (n = 236, ≤35 years) and control participants (n = 200) was carried out in the ethnic Amhara of rural North-West Ethiopia. We assessed their demographic and socioeconomic characteristics, and measured non-fasting C-peptide, diabetes-associated autoantibodies and HLA-DRB1 alleles. Leveraging genome-wide genotyping, we performed both a principal component analysis and, given the relatively modest sample size, a provisional genome-wide association study. Type 1 diabetes genetic risk scores were calculated to compare their genetic background with known European type 1 diabetes determinants. RESULTS: Patients presented with stunted growth and low BMI, and were insulin sensitive; only 15.3% had diabetes onset at ≤15 years. C-peptide levels were low but not absent. With clinical diabetes onset at ≤15, 16-25 and 26-35 years, 86.1%, 59.7% and 50.0% were autoantibody positive, respectively. Most had autoantibodies to GAD (GADA) as a single antibody; the prevalence of positivity for autoantibodies to IA-2 (IA-2A) and ZnT8 (ZnT8A) was low in all age groups. Principal component analysis showed that the Amhara genomes were distinct from modern European and other African genomes. HLA-DRB1*03:01 (p = 0.0014) and HLA-DRB1*04 (p = 0.0001) were positively associated with this form of diabetes, while HLA-DRB1*15 was protective (p < 0.0001). The mean type 1 diabetes genetic risk score (derived from European data) was higher in patients than control participants (p = 1.60 × 10-7). Interestingly, despite the modest sample size, autoantibody-positive patients revealed evidence of association with SNPs in the well-characterised MHC region, already known to explain half of type 1 diabetes heritability in Europeans. CONCLUSIONS/INTERPRETATION: The majority of patients with insulin-dependent diabetes in rural North-West Ethiopia have the immunogenetic characteristics of autoimmune type 1 diabetes. Phenotypic differences between type 1 diabetes in rural North-West Ethiopia and the industrialised world remain unexplained.

IL-1 Transcriptional Responses to Lipopolysaccharides Are Regulated by a Complex of RNA Binding Proteins.

The IL1A and IL1B genes lie in close proximity on chromosome 2 near the gene for their natural inhibitor, IL1RN Despite diverse functions, they are all three inducible through TLR4 signaling but with distinct kinetics. This study analyzed transcriptional induction kinetics, chromosome looping, and enhancer RNA production to understand the distinct regulation of these three genes in human cells. IL1A, IL1B, and IL1RN were rapidly induced after stimulation with LPS; however, IL1B mRNA production was less inhibitable by iBET151, suggesting it does not use pause-release regulation. Surprisingly, chromatin looping contacts between IL1A and IL1B were highly intermingled, although those of IL1RN were distinct, and we focused on comparing IL1A and IL1B transcriptional pathways. Our studies demonstrated that enhancer RNAs were produced from a subset of the regulatory regions, that they were critical for production of the mRNAs, and that they bound a diverse array of RNA binding proteins, including p300 but not CBP. We, furthermore, demonstrated that recruitment of p300 was dependent on MAPKs. Integrator is another RNA binding protein recruited to the promoters and enhancers, and its recruitment was more dependent on NF-κB than MAPKs. We found that integrator and NELF, an RNA polymerase II pausing protein, were associated with RNA in a manner that facilitated interaction. We conclude that IL1A and IL1B share many regulatory contacts, signaling pathways, and interactions with enhancer RNAs. A complex of protein interactions with enhancer RNAs emphasize the role of enhancer RNAs and the overall structural aspects of transcriptional regulation.

Genetic Discrimination Between LADA and Childhood-Onset Type 1 Diabetes Within the MHC.

OBJECTIVE: The MHC region harbors the strongest loci for latent autoimmune diabetes in adults (LADA); however, the strength of association is likely attenuated compared with that for childhood-onset type 1 diabetes. In this study, we recapitulate independent effects in the MHC class I region in a population with type 1 diabetes and then determine whether such conditioning in LADA yields potential genetic discriminators between the two subtypes within this region. RESEARCH DESIGN AND METHODS: Chromosome 6 was imputed using SNP2HLA, with conditional analysis performed in type 1 diabetes case subjects (n = 1,985) and control subjects (n = 2,219). The same approach was applied to a LADA cohort (n = 1,428) using population-based control subjects (n = 2,850) and in a separate replication cohort (656 type 1 diabetes case, 823 LADA case, and 3,218 control subjects). RESULTS: The strongest associations in the MHC class II region (rs3957146, ß [SE] = 1.44 [0.05]), as well as the independent effect of MHC class I genes, on type 1 diabetes risk, particularly HLA-B*39 (ß [SE] = 1.36 [0.17]), were confirmed. The conditional analysis in LADA versus control subjects showed significant association in the MHC class II region (rs3957146, ß [SE] = 1.14 [0.06]); however, we did not observe significant independent effects of MHC class I alleles in LADA. CONCLUSIONS: In LADA, the independent effects of MHC class I observed in type 1 diabetes were not observed after conditioning on the leading MHC class II associations, suggesting that the MHC class I association may be a genetic discriminator between LADA and childhood-onset type 1 diabetes.

Genetic and Epigenetic Fine Mapping of Complex Trait Associated Loci in the Human Liver.

Deciphering the impact of genetic variation on gene regulation is fundamental to understanding common, complex human diseases. Although histone modifications are important markers of gene regulatory elements of the genome, any specific histone modification has not been assayed in more than a few individuals in the human liver. As a result, the effects of genetic variation on histone modification states in the liver are poorly understood. Here, we generate the most comprehensive genome-wide dataset of two epigenetic marks, H3K4me3 and H3K27ac, and annotate thousands of putative regulatory elements in the human liver. We integrate these findings with genome-wide gene expression data collected from the same human liver tissues and high-resolution promoter-focused chromatin interaction maps collected from human liver-derived HepG2 cells. We demonstrate widespread functional consequences of natural genetic variation on putative regulatory element activity and gene expression levels. Leveraging these extensive datasets, we fine-map a total of 74 GWAS loci that have been associated with at least one complex phenotype. Our results reveal a repertoire of genes and regulatory mechanisms governing complex disease development and further the basic understanding of genetic and epigenetic regulation of gene expression in the human liver tissue.

Genome-scale Capture C promoter interactions implicate effector genes at GWAS loci for bone mineral density.

Osteoporosis is a devastating disease with an essential genetic component. GWAS have discovered genetic signals robustly associated with bone mineral density (BMD), but not the precise localization of effector genes. Here, we carry out physical and direct variant to gene mapping in human mesenchymal progenitor cell-derived osteoblasts employing a massively parallel, high resolution Capture C based method in order to simultaneously characterize the genome-wide interactions of all human promoters. By intersecting our Capture C and ATAC-seq data, we observe consistent contacts between candidate causal variants and putative target gene promoters in open chromatin for ~ 17% of the 273 BMD loci investigated. Knockdown of two novel implicated genes, ING3 at 'CPED1-WNT16' and EPDR1 at 'STARD3NL', inhibits osteoblastogenesis, while promoting adipogenesis. This approach therefore aids target discovery in osteoporosis, here on the example of two relevant genes involved in the fate determination of mesenchymal progenitors, and can be applied to other common genetic diseases.

First Genome-Wide Association Study of Latent Autoimmune Diabetes in Adults Reveals Novel Insights Linking Immune and Metabolic Diabetes.

OBJECTIVE: Latent autoimmune diabetes in adults (LADA) shares clinical features with both type 1 and type 2 diabetes; however, there is ongoing debate regarding the precise definition of LADA. Understanding its genetic basis is one potential strategy to gain insight into appropriate classification of this diabetes subtype. RESEARCH DESIGN AND METHODS: We performed the first genome-wide association study of LADA in case subjects of European ancestry versus population control subjects (n = 2,634 vs. 5,947) and compared against both case subjects with type 1 diabetes (n = 2,454 vs. 968) and type 2 diabetes (n = 2,779 vs. 10,396). RESULTS: The leading genetic signals were principally shared with type 1 diabetes, although we observed positive genetic correlations genome-wide with both type 1 and type 2 diabetes. Additionally, we observed a novel independent signal at the known type 1 diabetes locus harboring PFKFB3, encoding a regulator of glycolysis and insulin signaling in type 2 diabetes and inflammation and autophagy in autoimmune disease, as well as an attenuation of key type 1-associated HLA haplotype frequencies in LADA, suggesting that these are factors that distinguish childhood-onset type 1 diabetes from adult autoimmune diabetes. CONCLUSIONS: Our results support the need for further investigations of the genetic factors that distinguish forms of autoimmune diabetes as well as more precise classification strategies.

A Global Perspective of Latent Autoimmune Diabetes in Adults.

Latent autoimmune diabetes in adults (LADA) is characterized by the presence of islet autoantibodies and initial insulin independence, which can lead to misdiagnosis of type 2 diabetes (T2D). As such, understanding the genetic etiology of LADA could aid in more accurate diagnosis. However, there is ongoing debate regarding the exact definition of LADA, so understanding its impact in different populations when contrasted with type 1 diabetes (T1D) and T2D is one potential strategy to gain insight into its etiology. Unfortunately, the lack of consistent and thorough autoantibody screening around the world has hampered well-powered genetic studies of LADA. This review highlights recent genetic and epidemiological studies of LADA in diverse populations as well as the importance of autoantibody screening in facilitating future research.

Relative contribution of type 1 and type 2 diabetes loci to the genetic etiology of adult-onset, non-insulin-requiring autoimmune diabetes.

BACKGROUND: In adulthood, autoimmune diabetes can present as non-insulin-requiring diabetes, termed as 'latent autoimmune diabetes in adults' (LADA). In this study, we investigated established type 1 diabetes (T1D) and type 2 diabetes (T2D) genetic loci in a large cohort of LADA cases to assess where LADA is situated relative to these two well-characterized, classic forms of diabetes. METHODS: We tested the association of T1D and T2D GWAS-implicated loci in 978 LADA cases and 1057 non-diabetic controls of European ancestry using a linear mixed model. We then compared the associations of T1D and T2D loci between LADA and T1D and T2D cases, respectively. We quantified the difference in genetic risk between each given disease at each locus, and also calculated genetic risk scores to quantify how genetic liability to T1D and T2D distinguished LADA cases from controls. RESULTS: Overall, our results showed that LADA is genetically more similar to T1D, with the exception of an association at the T2D HNF1A locus. Several T1D loci were associated with LADA, including the major histocompatibility complex region, as well as at PTPN22, SH2B3, and INS. Contrary to previous studies, the key T2D risk allele at TCF7L2 (rs7903146-T) had a significantly lower frequency in LADA cases, suggesting that this locus does not play a role in LADA etiology. When constrained on antibody status, the similarity between LADA and T1D became more apparent; however, the HNF1A and TCF7L2 observations persisted. CONCLUSION: LADA is genetically closer to T1D than T2D, although the genetic load of T1D risk alleles is less than childhood-onset T1D, particularly at the major histocompatibility complex region, potentially accounting for the later disease onset. Our results show that the genetic spectrum of T1D extends into adult-onset diabetes, where it can clinically masquerade as T2D. Furthermore, T2D genetic risk plays a small role in LADA, with a degree of evidence for the HNF1A locus, highlighting the potential for genetic risk scores to contribute towards defining diabetes subtypes.